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1 Introduction 

Nonsymmetric elliptic eigenvalue problems have important physical background, 
such as convection-diffusion in fluid mechanics, environmental problems and so 
on. Thus, finite element methods for solving this problem become an important 
topic which has attracted the attention of mathematical and physical fields: [2] 
discussed a priori error estimates, u nn ® mi m m a posteriori error esti¬ 
mates and adaptive algorithms, [18] function value recovery algorithms, 
two level algorithms, nmi)] extrapolation methods, [5] an adaptive homotopy 
approach, etc. This paper turns to discuss finite element multilevel discretiza¬ 
tion based on local defect-correction. 

For elliptic boundary value problem, Xu and Zhou [21] combined two-grid 
finite element discretization scheme with the local defect-correction technique 
to propose a general and powerful parallel-computing technique. This technique 
has been used and developed by many scholars, for instance, successfully ap¬ 
plied to Stokes equation (see H3 mo, especially, Xu and Zhou (23], Dai and 
Zhou [8], Bi and Yang etc [4] developed this method and established local and 
parallel three-level finite element discretizations for symmetric elliptic singular 
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eigenvalue problems (including the electronic structure problems). 

In this paper, we further apply local defect-correction technique proposed by 
Xu and Zhou to nonsymmetric elliptic singular eigenvalue problems, our work 
has the following features. (1) We extend local and parallel three-level finite 
element discretizations for symmetric eigenvalue problems established by Dai 
and Zhou [8] to nonsymmetric eigenvalue problems. (2) Based on [4], we es¬ 
tablishes new multilevel finite element discretization by local refinement, this 
scheme repeatedly makes defect correction on finer and finer local meshes to 
make up for abrupt changes of local mesh size caused by three level scheme. 
And theoretical analysis and numerical experiments show that our schemes are 
simple and easy to carry out, and can be used to solve singular nonsymmet¬ 
ric eigenvalue problems. Numerical experiments show that, compared with the 
adaptive homotopy approach in [5], our algorithm seems to be more efficient. 
(3) For the nonsymmetric problems, based on the work of [201 El], we discuss 
the local error estimates of finite element approximations; its a new feature 
here that the estimates apply to the local domains containing corner points, see 
Lemmas 2.3-2.4 and Remark 2.2 in this paper. 

In this paper, regarding the basic theory of finite elements, we refer to 

ID El El US]. 


2 preliminaries 

Consider the nonsymmetric elliptic differential operator eigenvalue problem: 
d d 

Lu = — dj{aij{x)diU ) + ^bj ( x)djU + c(x)u = \m(x)u , in 12, (2.1) 

i,j =1 i =1 

u = 0, on dfl, (2.2) 

where 12 C R d , d > 2, be a polyhedral bounded domain with boundary dfl, 
® iU = * = " ' i d- 


Let 


a(u, v ) 
b(u, v) 



d 

+ bidiuv + cuv) dec, 

i=l 


n 

The variational form associated with m-m is given by: find A € C, 
u G #o(S2), ||u|| 0 = 1, satisfying 

a(u,v) = Xb(u,v), V v G Hq(Q). (2.3) 


Assume that a,ij,bi G IFi j00 (fI), c G Loo(I2) are given real or complex func¬ 
tions on Q, m £ Loo(12) is a given real function which is bounded below by a 
positive constant on 12. L is assumed to be uniformly strongly elliptic in 12, i.e., 
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there is a positive constant ao such that 
d d 

Re ^2 <k,j€i€j > a o^2 Va; e V(Ci,6, • ■ • ,£d) e (2.4) 

£,J = 1 2—1 

Let b = max |4(x)|. we assume without loss that i?e c > ao/2 + fr/(2ao) 
l<i<d;x€f2 

since adding a constant xm(x) to c{x) only shifts the eigenvalues. Under above 
assumptions, we have 

Re a(u,u) >^clq\\u\\\, Vu £ iL 1 (U); (2.5) 

and there are constants M\ and Mi such that 

|a(u,i>)| < Mi||u||i||i>||i, Mu, v G H 1 ^), (2.6) 

\b{u,v)\ < M 2 \\u\\ 0 \\v\\ 0 , \/u,v G L 2 (U). (2.7) 

For DcHoCll,we use D CC Slo to mean that dist(dD\dfl, <912o\<912) > 0. 
Assume that 7T/j( 12) = {r} is a mesh of f 1 with mesh-size function h(x) 
whose value is the diameter h T of the element r containing x, and h(Q) = 
rna xh(x) is the mesh diameter of 7^(12). We write h(Jl) as h for simplicity. 

Let 14(12) C G(12), defined on 7t/j(12), be a space of piecewise polynomials, and 
V°(f2) = T4(U) fl Hq (12). Given G C 11, we define nh(G) and 14(G) to be the 
restriction of 7pi(12) and 14(12) to G, respectively, and 

14°(G) = 14(G) n Hq(G), Vq(G) = {ve v°(n) ■. su PP « cc G}. 

For any G C 12 mentioned in this paper, we assume that it aligns with 7T/j. (12) 
when necessary. 

In this paper, G denotes a positive constant independent of h, which may 
not be the same constant in different places. For simplicity, we use the symbol 
cl < b to mean that a < Cb. 

We adopt the following assumptions in m for meshes and finite element 
space. 

(AO) There exists v > 1 such that h(fl) 1 ' < h(x), Vx G 12. 

(Al) There exists r > 1 such that for w G H^fl) G iJ 1+t (12), 

inf (\\h~ 1 (w - u)|| 0 + ||w - n||i) < /i‘||u;||i +t , 0 < t < r. 

v&V°(Cl) 

(A2) Inverse Estimate. For any v G 14(12 0 ), ||r>||i,n 0 ^ ||^ _1 ^||o,n 0 ■ 

(A3) Superapproximation. For G C 12o, let w G G°°(12) with supp cc CC G, 
then for any w G 14(G), ic|dGndn = 0, there exists v G V 0 h (G) such that 

\\h~ l (ww - u)||i,g < IM]i,g- 

Let 7r^(12) consist of shape-regular simplices and (AO) hold, and let 14 (12) C 
G(12) be a space of piecewise polynomials of degree < r defined on 7 t/i( 12), then 
from {21] we know that (Al)-(A3) are valid for this 14(12). 
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The finite element approximation of m is given by: find A h £ C, Uh £ 
V°(D), Iloilo = 1, satisfying 

a{u h ,v) = \ h b(u h ,v), Vr € V ft °(0). (2.8) 

Thanks to [T], we know the adjoint problem of (l2.lH - fl2.2D is: 
d d 

L*u* = — ^ di(aijdjU*) — ^ di(biU*) + cu* = \*mu *, infl, (2.9) 

i,j=l i= 1 


u* = 0, ondfl. (210) 

The corresponding variational form and discrete variational form of (I2.9D - (12. 10D 
are given by: find A * £ C, u* £ ||w*||o = 1, satisfying 

a(u,u*) =Yb(v,u*), V v £ Hq (12); (2.11) 

find \* h £ C,u* h £ 14, K||o = 1, satisfying 

a{v,u* h ) = >* h b(v,u* h ), Vv£V°(n). (2.12) 

Note that the primal and dual eigenvalues are connected via A = A* and 
Aft = K, 

Throughout this paper, we will assume that m-m hold. Thus from 


Lax-Milgram theorem we know the source problem associated with (12.311 and 
(12.1111 admits an unique solution, respectively. The discrete source problem as¬ 
sociated (12.81) and (12 121) admits an unique solution, respectively. 

Define the solution operator T : L 2 (£l) —> Hq(£1) and Th : L 2 (£l) —> V)°(0) 
as follows: 


a{Tg,v) = b(g,v), \/v £ (2.13) 

a(T h g,v) = b(g,v), Vu £ V°(£l). (2.14) 

And (12.311 and (12.811 have the equivalent operator form (12.1511 and (12.1611 . 
respectively. 

Tu = \~ 1 u, (2.15) 

ThUh = A jj^Uft- (2-16) 

Define the solution operator T* : L 2 (Q) — v Hq{£1) and T£ : L 2 (Cl) —> V®{VL) 
satisfying 

a(v,T*f) = b(v, /), Vue (2.17) 

a(v,T*f) = b(v,f), Vue 14°(D). (2.18) 

And (12.111) and (12.1211 have the equivalent operator forms (12.191) and (12.201) . 
respectively. 

T*u* = A*"V, (2.19) 

T* h u* h = A r 1 ^. (2.20) 
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It can be proved that T is completely continuous, and T* is the adjoint 
operator of T in the sense of inner product b(-, ■). In fact, 

b(Tu,v ) = a(Tu,T*v) = b(u,T*v), Vu,v G P 2 (fl), 
b(T h u, v ) = a(T h u, T£v) = b(u, T£v), Vm, v G L 2 (VL). 

We need the following regularity assumption. For any / G L 2 (VL), Tf G 
Hq(£1) PI H 1+11 (fl) and T* f G H^fl) n H 1+l2 (fl) satisfying 

||T/|| 1+71 <Co||/||o, (2.21) 

||T*/|| 1+72 < Ca||/||o. (2.22) 

According to m and the section 5.5 in [3J, the above assumption is reasonable. 

For some G C SI, we need the following local regularity assumption. 

R(G). For any / G L 2 (G), there exists a <j) G Pq(G) H H 1+1i (G) satisfying 

v) = b(f, v), Vv G Hl(G), 

and 


||<(>||i+7i,g < CgII/II-i+ti.g- (2.23) 

For any g G L 2 (G), there exists ai^G Hq (G) n H l+l2 (G) satisfying 
a{v,ip) = b(v,g), Vu G Hq(G), 


and 


IMIi+72,g < C' G ||g||_i +72jG . (2.24) 

Where Cn, C G are two priori constants, and not necessarily the same at differ¬ 
ent places. 

Define the Ritz projection Ph : Hq(£1) —>• and : Hq(£ 1) V®(Cl) 

by 

a{u — PhU, v) = 0, and a(v, u — P^n) = 0, G V® (12). (2.25) 

Then T,, = P^T, T* = P*T* (see [I]). 

Let M (A) be the space spanned by all generalized eigenfunctions correspond¬ 
ing to A of T, Mh{ A) be the space spanned by all generalized eigenfunctions 
corresponding to all eigenvalues of Th that converge to A. In view of the adjoint 
problem (12.11[) and (12.121) . the definitions of M*( A*) and M£(A*) are analogous 
to M( A) and M/,(A). 

In this paper, we suppose that A is an eigenvalue of (12.31) with the algebraic 
multiplicity q and the ascent a = 1. Then A* = A be eigenvalue of (12.111) . M( A) 
and M*( A*) are all eigenfunction space. 

Let Ah be the eigenvalue of (12.81) which converges to A, let A£ = Ah, and 
M*( A£) be the generalized eigenfunction space corresponding to the eigenvalue 
A*ofT*. 
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Remark 2.1. Obviously, it’s difficult to determine the ascent a of the 
eigenvalue A of (12.3D theoretically. But one could easily find that when the 
ascents of the eigenvalues of (12.81) . which converge to the same eigenvalue A 
of (12.31) . are all equal to 1, one can conclude that the ascent a = 1 from the 
standard theory of spectral approximation. And the ascents of eigenvalues of 
(12.81) can be determined by computation. 

We also need the lemma as follows (see [I6lf25l ): 

Lemma 2.1. Let (A, u) be an eigenpair of (12.31) . and (A* = A, u*) be the 
associated eigenpair of the adjoint problem (12.111) . Then for all w, w* £ Hq (12), 
b(w , w*) / 0, 

a{w,w*) _ _ a(w u, w* u*) _ b(w -u,w* - it*) 

b(w,w*) b(w,w*) b(w,w*) 

Proof, see [16J [25]. □ 

The a priori error estimates of the finite element approximations (12.81) and 
(12.121) can refer to pQ [2]. 

Lemma 2.2. Assume that M(A) C H r+S (Q), M*( A*) C H r+S2 (Q) (0 < 
s, S 2 < 1). Then 

|A/i — A| < h r+s ~ 1+r+S2 ~ 1 - 1 (2.27) 

let Uh £ Mh( A) with Hu^Ho = 1, then there is u G M(A) such that 

\\u h -u\\ 1 <h r+ -\ (2.28) 

\\uh-u\\o < h r+s ~ 1+ry2 - (2.29) 

let G M^{ A*) with || 0 = 1, then there is u* G M*( A*) such that 

K-«li^* r+ ' a " 1 ; (2-30) 

K-«lo<* r+ ' a_1+71 . (2.31) 


Proof, see [Tj. □ 

etc. studied the local behavior of finite element. The following 
Lemma 2.3 is a simple generalization of Lemma 3.2 in [21] . We can easily prove 
this Lemma by the same argument as that of Lemma 3.2 in [2lj . 

Lemma 2.3. Suppose that / G H~ l (Cl) and G CC fto C 12. If w G I4(flo), 
w|annan 0 = 0> satisfies 

a(w, v) = f(v), Vv G Vq(Qo), (2.32) 

then 

IMIi.g $ lk||o,o 0 + ||/||-i,n 0 j (2.33) 

where 


11/11—1,^0 = SU P fW- 

0efl'j(no),ll0l|i,ri o =i 
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Proof. Let p > 2v — 1 be an integer, and let 

D CC ttp CC fip-1 CC---CC fil CC n 0 . 

Choose Di C U satisfying D CC D± CC fl p and w G C°°(Cl) such that 
supp w CC ilj, and w = 1 on D\. Then, from (A3), there exists v G V 0 h (fl p ) 
such that 


||w 2 w - u||i,f2 p < hn 0 |klli.n p , 


so we have 


a{w, uj 2 w -v)< hn 0 |klli,n p 

and 


(2.34) 


\f(v)\ < ll/ll-i,n 0 IMIi,n p 

< ||/||-i,n 0 (ftfi 0 |M|i,fi p + IkHli.n)- (2.35) 

Since v G V 0 h ({l p ) C V 0 h (fio), (12.321) implies 

a(w,uj 2 w) = a(w,co 2 w — v) + f(v). (2.36) 

d 

Let ao (u, v) = f dijdiudjV. We can be derived from proof of Lemma 3.1 in 

i,j= 1 

[2l] that if Q 0 C C R d (d = 2,3), w G C'°°(fi) , supp ui CC flo, then 

ao(cow, unv) < a(w, w 2 u>) + ]kllo,n 0 ’ Vw G Hq(U). (2.37) 

It follows from (I2.34I) - (I2.37I) that 

IkHli.n ^ a 0 (uw,u. >w) < ci(w,uj 2 w) + |kllo.n 0 
= a(w,uj 2 w -v) + |kllo,n 0 + / (v) 

^fiolklli.fip + lkllo.no + ll/ll-i,n 0 (hn 0 ||u;||i,n p + ||aJto||i,n), 

thus 

Iklli.D < h 1 ^ |M|i,n p + Ikllo.fio + ||/l|-i,n 0 - (2.38) 

Similarly, we can get 

IHIi.fii - + lkllo.no + ll/ll-i.Oo, 3 = 1,2,-•• ,p. (2.39) 

By using (12.381) and (12.391) . we get from (AO) and (A2) that 

Iklli.D < < +1)/2 |klli,n 0 + lkllo.no + ||/|| —i,n 0 

< < +1)/2 |k- 1 ^|| 0> n o + lkllo.no + ll/ll—i,n 0 

< |kllo.n 0 + ll/ll—i,f 2 0 - 

This completes the proof. □ 
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Lemma 2.4. Suppose that G CC flo C SI. Then the following estimates 
are valid: 

h? 2 \\u - P h u\\i q + ||u - Ph,u\\ 0 'n < h 12 inf ||u - d|i,n, (2.40) 

\\u - Ph,u\\i,G < inf ||u - u||iA> + h l2 \\u - P h u\\i ,n- (2.41) 

vev°(n) 

Proof. For proof of (12.401) cf. [3, 7j, for proof of (12.411) cf. Theorem 3.4 in 

m- □ 


Remark 2.2. In |21| . the condition Superapproximation is given as follows. 
A. 3. Superapproximation. For G C Sl 0 , let w G Co°(f2) with supp u> CC G. 
Then for any w € 14(G), there exists v G V 0 h (G) such that 

II h^iuw - u)||i jG < ||uj||i,g- 

In the proof of Lemma 3.2 in m the authors choose D\ C satisfying D CC 
D\ CC Sip and ui G Cq°(SI) such that w = 1 on D\ and supp u ; CC Sl p . 

This paper just makes a minor modification, so that the theory of the local error 
estimates built in [21] applies to the local domains containing the corner points, 
see Lemma 2.3 and Lemma 2.4. 


3 Multilevel discretizations based on local defect- 
correction 

Consider the eigenvalue problem (12.31) which has an isolated singular point z G 
SI (e.g., see Figure 3.1). 

Let D CC SI be a given subdomain containing the singular point z, and we 
introduce domains 

SI D Hi D Sl 2 D • • • D Sli DO D. 

Let tth{SI) be a shape-regular grid, which is made up of simplices, with 
size H G (0,1), 7r w (Sl) be a refined mesoscopic shape-regular grid(from ttb(SI)) 
and TTh(Sli) be a locally refined grid (from 7T/ Z . i _ 1 (r2^_i)) that satisfy h_i = H , 
h o = w, hi -C hi -1 (i = 0,1, ■ • * ,1). (Figure 3.1 shows tvh{SI), tt w (SI) and 
7r^(f2i)). Let Vjj(Sl), V®(Sl), and {V®.(Sli)}[ be finite element spaces of degree 
less than or equal to r defined on tth{SI), 7r w (Sl) and {7r/ li (f2*)}i) 5 respectively. 



Figure 3.1 





























































Based on Algorithm Bq in [5] we establish the following three-level discretiza¬ 
tion scheme. 

Scheme 3.1 (Three-level discretizations based on local defect-correction.). 
Step 1. Solve (12.31) on a globally coarse grid 7Th(12): find A h G C. uh G 
Vff(fi) such that ||w_f/||o = 1 and 

a(u H , v) = \ H b(u H , v), Vu G 


Let A^ = A h, and find u * H G with ||m^||o = 1 such that \b{uH,u* H )\ 

has a positive lower bound uniformly with respect to H (see section 5.1). 

Step 2. Solve two linear boundary value problems on a globally mesoscopic 
grid n w (Q): find u w G V!^(fi) such that 

a(u u \v) = X H b(u H ,v), Vu G 1^(12); 


Hnd u w * G Ku(n) such that 

a(v,u w *) = \ H b(v,u* H ), VuG 1^(12). 

Then compute the Rayleigh quotient A u ’ = ■ 

Step 3. Solve two linear boundary value problems on a locally fine grid 
tt/h^i): find e hl G V^Qi) such that 


find e hl * 


Step 


a(e h \v) = X w b{u w ,v)-a{u w 

», VuG <(12,); 

(3.1) 

G VJ^fii) such that 



a(v, e hl *) = A w b(v, u w *) — a(v, 

u w % VuG <(120. 

(3.2) 

4. Set 



_{ u w + e ftl 

U ~ \ u w 

on 12i, 
m 12 \ 12i 

(3.3) 

n , w , hl * f u w * + e h '* 

U ~ \ u w * 

on f2i, 
in 12 \ 12i 

(3.4) 


And compute the Rayleigh quotient 


^w,h i 


a{u w ’ hl ,u w ’ hl *) 
b(u w ’ hl ,u w ’ h i*) ’ 


yw,h x * _ y i w,h 1 


(3.5) 


We use (X w ’ hl ,u w,hl ) and (X w,hl *,u w,hl *) obtained by Scheme 3.1 as the 
approximate eigenpair of (12.31) and (12.111) . respectively. 

It is obvious that (A w ,u w ) and (A w ,u w *) in Scheme 3.1 can be viewed as ap¬ 
proximate eigenpairs obtained by the two-grid discretization scheme in PS [25] 
from 7r/f(Sl) and n w (fl). 


Using Scheme 3.1, abrupt changes of mesh size can appear near <912 1 . Influ¬ 
enced by the technique on the transition layer proposed by |U, we repeatedly 
use the local defect-correction technique to establish the following multilevel 
discretization scheme. 

Scheme 3.2 (multilevel discretizations based on local defect-correction.). 
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Step 1 .The same as that of Step 1 of Scheme 3.1. 

Step 2 .The same as that of Step 2 of Scheme 3.1. 

Step 3. u w ’ h ° <= u w , X w ’ h ° <= X w , u w ’ ho * <= u w *, X w ’ h °* <s= X w *. 

Step 4. For i = 1,2, • • • ,1, execute Step 5 and Step 6. 

Step 5. Solve linear boundary value problems on locally fine grid 7^(12*): 
find e hi G V®. (12*) such that 

a(e h \v) = X w ’ hi - 1 b{u u '’ hi -\v)-a(u w ’ hi - 1 ,v), Vv G V£(fii); (3.6) 
find e hi * G such that 

a{v, e hi *) = X w ’ hi ~ 1 *b(v, u w ’ hi ~ 1 *) - a(v, u w ’ h Vu G (12*)- (3.7) 


Step 6. Set 


u 


w,hi 


,hi -1 
,hi-i 


on 12*, 
in 12 \ 12* 


(3.8) 


And compute 


w,hi 


u wM-i* + e hi* 

yWjki—l* 


on 12*, 
in 12 \ 12* 


X w ’ hi = 


a(u w,hi , 

b(u w ’ h \ 


u w ’ hi *) 
u w < h <*)’ 


yv,hi* — \w,hi 


(3.9) 


(3.10) 


We use (X w,hl ,u w,h ‘) and (X w,hl *, u w,hl *) obtained by Scheme 3.2 as the 
approximate eigenpair of (12.311 and (12.111) . respectively. 


4 Theoretical Analysis 

Next we shall discuss the error estimates of Scheme 3.1 and Scheme 3.2. 

In our analysis, we introduce an auxiliary grid 717 ** (12) which is defined 
globally, and denote the piecewise polynomials space of degree < r by V® ( 12 ) 
[i = 1, 2, ■ • • , V). We also assume that (12*) and V°. (12*) are the restrictions 
of 7 r/u(fi) and F°.(lI) to 12 *, respectively, and 

vg(fi) c v£(fi) c <(i2) c vft(fi) c • • • c v£(n). 

For D and 12,; stated at the beginning of section 3, let Gi C 12 and F C 12 
satisfy D CC F CC Gi CC 12* (i = 1,2,--- , 2). 

Theorem 4.1. Assume that M(A) C 7/,) (12) (~l iF +s (12) fl iF +1 (12 \ D) 
and (1 < r + s, 0 < s < 1), M*(X*) C n FF +S 2 (12) n £F +1 (12 \F) and 

(1 < r + S 2 , 0 < S 2 < 1), and H is properly small. Then there exists u G M(X) 


G M*{ A*) such that 

II u w - m||i 

< 

H r+s-l+ l2 +w r+ S -l j 

(4.1) 

\\u w ~ u|| 0 

< 

jjr+s— 1+72 

(4.2) 

IK -«lli,n\3? 

< 

jyr-M-l+T. +w r, 

(4.3) 

Hu™* -u*||i 

< 

JJH--2-1+7! +U ,r+ S 2 -1 ) 

(4.4) 

IK*-«*llo 

< 

j^r+s 2 - l+7i 

(4.5) 

\\ uW * ~ M *lli,n\P 

< 

H r+ S2 -l+ 7 i +w r, 

(4.6) 

| X w - X | 

< 

yj 2 r+s+s 2 - 2 + 7 i + 72 _|_ w 2r+s+s 2 -2 

(4.7) 
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Proof. Let u G M{ A) and u* G M*(X*) such that u — uh and u* —u* H both 
satisfy Lemma 2.2. From (12.131) . (12.141) . Step 2 of Scheme 3.1, (12.151) . Lemma 
2.2 and Lemma 2.4, we derive that 

II u w - u\h = \\\hT w u h - ATu||i 

< || A h T w u h - Ar,„u||i + HAT^m - ATw||i 

< IIA HUH - Au||o + A|| P w Tu - Tu\\i 

< jjr+s- 1+72 _|_ w r + s ~ 1 

then dUD follows. By Lemma 2.2 and Lemma 2.4, 

ll M — u IIi,q\f ~ IIA huh — Au||o + X\\P w Tu — TwII^Q^p 
< H r+S ~ 1+ 72 + w r , 
then (14.31) follows. By calculation, 

IK - u|| 0 = || A h T w u h - ATu|| 0 

< || A h T w u h - ATu,m || 0 + ||ATK - ATu|| 0 

< ||A H u H - Au|| 0 + A|| P w Tu - Tw|| 0 

< jjr+s- 1+72 _|_ w r + s ~ 1+72 

< }J rJrS ~ 1+72 

then (14.21) follows. 

Similarly we can prove (EES, (EES and (14.61) . From (12.261) . we have 

_ A = a(u w -u,vP*-u*) _ ^ b(u w — u, u w * — u*) 

b(u w ,u w *) b(u w ,u w *) ' 1 ' ’ 

Note that uh and u w just approximate the same eigenfuntion u, u* H and u w * 
approximate the same adjoint eigenfuntion u*, \b(un,u* H )\ has a positive lower 
bound uniformly with respect to H , therefore b(u w ,u w *) has a positive lower 
bound uniformly. Combining gm) , q, dm, (EEs and m yields dm- □ 

The following Theorem 4.2 is a critical result in this paper, which develops 
the results of Theorem 3.3 in [8l. 

Theorem 4.2. Assume that holds (i = 1,2, • • • , l), u G M(A) and 

u* G M*(X*). Then 

I \u wM - P hl u\kn < ||u - P hl u || 0 ,n, + - u w,hl - 1 l|i,n, 

+ || Aw - X wM -2u wM -i\\ 0tQl + || /«-!„»,/.,_i _ Au || 0 

+ ||u 1 1 — -fh-i u lli,n\Gr 4” H M ’ 1 1 — u lli,f)i\F> ^ — (4-9) 


\\u wM * - P fc >||i,n < ||«* - P* hl n lo,n, + 

i || \ * \ w,hi_o* „w,hi_ o* II i || \ w,hi_-y *„w,hi_ i * \„*m 

+ ||A U - A ’ 1 2 U ’ 1 2 ||o,ft, + ||A ’ 1 1 U ’ 1 1 — AU ||o 


U 


w,hi_i* 


~ ^« u *lli,f2\G7 + \\u wM - 1 * - u*||i A \y, l > 1.(4.10) 
Proof. Due to the inequality 

I \u wM - ^.ulli.n < \W wM - P hl u\\i,D + ^ PhM\i, Gi \d 

+ \\u w ’ h ‘-P hl u\\ lt n\^ (4-H) 
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we shall estimate \\u w ’ hl -P hl u\\ i >D , \\u w,hl -Ph l u\\i,G,\D’ and W u ' w,hl ~ p hi u \\i,n\G l 
, respectively. 

First, we proceed to estimate \\u w ’ hl — Ph,u\\i t D- From (13.81) . (13.61) and (12.251) 
we derive 

a(u w,h 1 — P hl u, v) = a(u w,h ‘,v) — a(Ph t u, v ) = a{u w ’ hl ~ l + e hl ,v) — a(u, v ) 

= X w ’ h ‘- 1 b(u w ’ h, - 1 ,v)-Xb(u,v), VvG V£(fii)- (4.12) 

It is obvious that 

X w ’ h, - 1 b(u w ’ hl - 1 ,v) - Xb(u,v) 

= (X w ’ hl ~ 1 - X)b(u, v) + X w,hl - 1 b(u u '’ hl - 1 - u, v), Vv G ^(n)j(4.13) 
which together with (14.121) yields 

a( u w ' h ‘ - P hl u,v) = (X w ’ h ‘~ 1 - X)b(u,v) + X w ’ h ‘- 1 b(u w ’ h ‘~ 1 -u,v), Vv G Vg^Qi). 

Since Lg (fl;) C V® (i\) , thus, from the above formula and Lemma 2.3 we 
deduce that 

|| U w ’ h ‘ - P hl u\\i,D < \\u w ’ hl - P hl u\\ o.o, + \X w ’ h, ~ 1 - X\ + || u w ’ h ‘ - u||oa(4.14) 

By calculation, we have 

IK’"* -PnAWn, < Wu^- 1 -PhMkn, + He^lloA 

< ||m — Phiu\\o,sh + ||u — 1 || 0 ,n, + ||e hl ||o,n,, 

substituting the above relation into (14.141) we obtain 
\\ U W ^ - P hl u\\l,D 

< \\wM-i _ A| + - U || 0A + || U - i\,u||oA + lie"' 110,0,(4.15) 

To estimate \\e hl ||o,Oi, we use the Aubin-Nitsche duality argument. For any 
given / G L- 2 {p,i), consider the boundary value problem: find ip G such 

that 

a(v, ip) = b(v, f) Vv G Hq (Cli). (4.16) 

Let ip be the generalized solution of (14.161) . tph t and iph^, be finite element 
solutions of (14.161) in V^(f2/) and V£ I _ i (f2i), respectively. Then, 

Ib-^hJIiA ||oa, < /i?!!il|/||o,o r (4.17) 

From (13.61) and (13.81) we get 

a(u w ’ h ‘,ip hl ) = X w ’ hl - 1 b(u w,hl - 1 ,tp hl ), 
thus by the definitions of ip, iphi and e hl , we deduce that 

b(e hl , f) = a(e h ‘,ip) = a(e hl ,ip hl ) = a(u w ’ h ‘ - u w ’ h, ~ 1 , ip hl ) 

= a(P hl u - u w ' h ‘- 1 ,ip hl ) + a(u w ’ hl ,ip hl ) - a(P hl u,<p hl ) 

= a(P hl u - u w ' h ‘~\ip hl ) + X w ' h ‘-'b(u w ’ h ‘-',<p hl ) - A b(u, <p hl ) 

= a(P hl u - u w ’ h ‘- 1 ,<p hl -<p) + a(P h ,u - u w ’ h ‘- 1 ,tp- tp h ,_J 
+ a(P hl u - u w ' hl ~ 1 ,ip hl _ 1 ) + X w ’ hl - 1 b(u w ’ hl - 1 ,ip hl ) - Xb(u,ip hl ) 

< h^WPb'U - u w ^ || 1A 11/110,0, + a(P hl u 
+ X w ’ hl ->b(u w ’ h ‘~\ip hl ) - Xb(u,tp hl ). (4.18) 


12 












Step 2 of Scheme 3.2 shows that 

_( n .w.hn \ \w,h-\T( n .w,h-\ \ 

d\U i^Pho) — ^ 0 \U 5 ^Phoji 

namely, for 1 = 1, 

a{u w > h '~\ip hl - J = ^ w ’ hl ~ 2 K uW ' h, ~ 2 j l Phi- 1 ), 

for l > 1, the above formula follows from (13.611 and (13.81) . Therefore, 

a(P hl u - u w ’ h ‘- 1 ,ip hl _ 1 ) = a(u - u w ’ hl ~ 1 ,<p hl _ 1 ) 

= A6(u, <£/,,_,) - a{u w ’ hl - 1 ,<p hl _J 

= A6( U ,93ft,_ 1 )-A“'> ft ‘- 2 6(t i “'’ , “- 2 ^ft,_ 1 ) 

^IIAu-A^-^^-IIoaII/IIoa- 

It is clear that 

|A’"- fc «- 1 6(u”’> fc «-S ^ft.) - A6(w, ^ft,)| < \\\™’ h ‘-iu w ’ h ‘-' - A«||o,n,II/II oa- 
Substituting the above two formulae into (14.18[l . we derive 

I b(e h ‘,f)\ < (^rll^u-u^'-MliA + ||Aw - A™> h «-’«“- h '-’|| 0 , ni 

+ ||A^- 1 ^- 1 -Au|| 0 )||/||o,o i . 

Thus, we get 

lle^lloA < h^WP^u - u w ’ h ‘-i\\ i,o, + ||Aw - || 0 ,n, 

+ ||A’"- h «- 1 u“’. fc «-i - Aw|| 0 . (4.19) 

Substituting (14.191) into (14.151) . we obtain 

II n w ’ h ‘ - P hl u\h,D < ||u - Pft^HoA + h^JP^u - u w,hl ~ 1 ||i,n, 

+ ||Aw - X w ’ hl - 2 u w ’ hl - 2 \\o,ni + \X w ’ hl ~ 1 - A| + \\u w ’ hl ~' - w|| o .(4.20) 
Similarly, since (G; \ D) CC f h, we deduce 

\\u w ’ hl - P hl u\\i,D < ||u - Pft^HoA + h^JP^u - w^-IIia 

+ ||Aw - A“> , ‘*-=>u w > fc «-=>||o i n I + \\ w ’ hl ~ 1 - A| + \\u w ' hl ~' - w|| 0 .(4.21) 

The remainder is to analyze \\u w ’ hl — Ph l u\\ l From (13.81) . we see that 

h w ' hl - PhMknxu, = -Pft,w||i MTTi , 

which leads to 

ll 11 ' 1 — 

— \\u W ’ hl — -PftiW|| i n \Q, + ||w ’ 1 1 — Pft,w||i n ,\Qj-+ ||e -I || 1> f 2 ,\G 7 

< -Pft^H^Gl + ll^'lknAGf ( 4 - 22 ) 

It follows from (13.61) . (12.31) and (14.131) that 

a(e hl ,v) = \ w ’ h ‘- 1 b(u w,hl - 1 ,v) — a(u w ’ h ‘~ 1 ,v) — Ab(u,v ) + a(u,v) 

= (\ w ’ h i-i - X)b(u, v) + A w ’ hl - 1 b(u w ’ hl ~ 1 - w, v) 

— a(u w ’ hl ~ 1 — w, v), Vw G T4°(f h), 
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then, by Lemma 2.3, we have 


\ ehl Hi, n,\G, ~ 


o h l\ 


\e-‘ Ho,+ \ X - “HiA\ F’ ( 4 ' 23 ) 


where F C XI satisfies fiCCfCCGi. Substituting (14.231) into (14.221) we get 


i’ 0 -'* 1 - P, 




< 


- ^.«lli,n\sr + Iloilo,n A P + lA ^'- 1 - Al + 


u 


w,h i _ 1 


h,n,\F- 


It follows from substituting (14.191) into the above inequality that 

\\u wM - P»M\i, a \7T, % IK’'*'- 1 ^ P*M\i,n\(T, + W-rWP^u - 
+ ||Aw - A u '’ ft 'i- 2 u"' ,/i ' i - 2 || 0 i n i + ||A v>M-i u v>,h,-i _ a «|| 0 

+ \\^- 1 -X\ + \\u^ h ^-u\\ ltaA¥ . (4.24) 


Combining (14.241) . (14.201) . (14.211) and (14.111) . finally, we obtain (14.91) . 
We can prove (14.101) by using the similar argument. □ 


Theorem 4.3. Assume that the conditions of Theorem 4.1 hold. Then 
there exists u G M(A) and u* G M*{ A*) such that 

\\u w ’ hl - u||i, n < h r + s ~ l +w r + H r+S ~ 1+12 , 

|| u w ’ hl - u||o,n < w r + H r+S ~ 1+ ^, 

\\u w ’ hl ~ u lli,fi\F < w r + H r+S ~ 1+ ^ 2 , 

\\u w ’ hl * _ u *|| i n < h^ 82 - 1 + w r + H r+S2 ~ 1+ ^, 

\\u w ’ hl * - u*|| 0 ,fi <w r + H r+S2 ~ 1+ i i, 

(I u w ’ hl * - u*|| i n ^ < w r + H r+S2 ~ 1+7 \ 

_ ^| < ^2r+s+s 2 -2 yj2r _|_ ^y 2 r+s+s 2 - 2 + 7 i+ 7 2 

Proof. Let u G M{ A) and u* G M*(X*) such that u — uh and u* — u* H 
both satisfy Lemma 2.2. In Theorem 4.2, choose l = 1, = H, ho = w, 

u w,h 0 = u w^ X w,h 0 = X w, u w,h^ = Uh ^ X w,h -i = Xh , then we get 

|| U W ’ kl - PftiUHl.fi < ||U - PftiUHo.fi! + W 72 ||PftiU - U W || l.fii 
+ ||Au — A_fftUftf||o.fii + HA^u™ — Au||o 

+ ||u — P/ti u lli,fi\Gr H u — M Hi,fii\F- (4.32) 

Using Lemma 2.4, Theorem 4.1, Lemma 2.2 to estimate the terms at the right 
hand side of the above formula gives 

| \u W ’ hl - PftiUHi.fi < h\ +S - 1+12 + W l2 W r+S ~ 1 + LT+s-1+7 2 + w r+s- 1+ 72 

+ (uC +s - 1+72 + w r ) + (uP+ s - 1+72 + w r ) < H r+S ~ 1+ 1 2 + w r . (4.33) 


(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 


Combining (12.401) and (12.411) yields (14.251) . (14.261) and (14.271) . By the same 
argument we can prove (14.281) . (14.291) and (14.301) . From (12.261) . we have 


X w ’ h1 - A = 


a(u w,hl - u,u w ’ hl * -u*) 
b(u w ’ hl ,u w ’ hl *) 


-X 


b{u w ’ hl — u, u w,hl * -u*) 
b(u w ’ h 1 ,u w ’ hl *) 


(4.34) 


Note that uh and u w,h 1 just approximate the same eigenfunction u, u* H and u w ’ hl * 
approximate the same adjoint eigenfunction u*, \b(uH, u* H ) \ has a positive lower 
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bound uniformly with respect to H , thus b(u w ’ hl , u w,hl *) has a positive lower 
bound uniformly. Combining (14.251) . (14.261) . (14.281) . (14.291) and (14.341) yields (14.311) . □ 


For convenient argument, we assume S 2 = s, 7i = 72 = 7 in the following 
Theorem. 

Theorem 4.4. Under the conditions of Theorem 4.1, we further assume 
that holds (i = 1, 2, ■ • • , /), and 

w r = 0(H r+s ~ 1+7 ), h\ +s ~ l > H r+S ~ 1+ U (4.35) 

Then there exists u € M(A) and it* £ M*(A*) such that 

,w,hi 




hs-l 

w||o,n < ff r+s “ 1+7 , 






< hr+ 


s—1 




u w,hi*_ u * || on < 


- A| < h 2r+2s 


(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 


Proof. Let u £ M(A) and u* £ M*(A*), such that u — uh and w* — 
both satisfy Lemma 2.2. The proof of (14.36I) - (I4.42I) is completed by induction. 
When l = 1, Scheme 3.2 is actually Scheme 3.1. Hence, from Theorem 4.1, 
Theorem 4.3 and (14.351) we know that (14.36[) - (14.42[) hold for l = 0,1. 

Suppose (14.361) - (14.421) hold for l — 2, l — 1, i.e., 


and 



(4.43) 

\\u w ’ h ‘- 2 - u||o, n < H r+a ~ 1+ ~<, 

(4.44) 

^ u W,hl-2 _ U i < H r+S ~ X+1 , 

(4.45) 


(4.46) 

\\u w ’ hl ~ 2 * - u*|| 0 ,n < H r+S ~ 1+1 , 

(4.47) 

\\u w ’ hl ~ 2 * - u*|| 1M ^ < H r+S ~ 1+ 7 , 

(4.48) 

\\ w ’ hl ~ 2 - A| < h 2r _+ 2s ~ 2 - 

(4.49) 


(4.50) 

\\u w ' hl -' - u|| 0 , n < ^+*- 1 + 7 , 

(4.51) 


(4.52) 

II \u.hi i* „,*|| cr' ir+s-l 

||^ ^ 111,0 — l ? 

(4.53) 

\\u wM ~^*- u *\\ 0i n<H r+s - 1+r y, 

(4.54) 

|| w ^«- i *_ u *|| iM _<^+ S -i+ 7 , 

(4.55) 

lA^.U-i _ A | < h 2r+2s-2 

(4.56) 


Next we shall prove that (14.361) - (14.421) hold for l. Using the above formula and 
Lemma 2.4 to estimate the terms at the right hand side of (14.91) gives 

|| u w ’ h ‘ - P hl u\\i,n < h \ +s ~ 1+1 + h'l_ 1 (h r + s ~ l + h^- 1 ) + H r+S ~ 1+1 
+ H r+S - 1+1 + (H r+S - 1+ ^ + h\) + H r+S ~ 1+1 < H r+S ~ 1+7 (4.57) 
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Combining (12.401) . (12.411) and (14.571) yields (14.361) . (14.371) and (14.381) . By the 
same argument we can prove (14.391) . (14.401) and (14.411) . From (12.261) . we have 


X wM _ A = 




,w,hi _ , 




- «*) x K 


b(u w ’ hl ,u w ' h i*) 


-A- 


VI,hi _ , 


l w ' hl * - U* 


b(u w ’ h ‘,u w ’ h '*) 


(4.58) 


Using the similar argument as that of Theorem 4.3 we know that b(u w ’ h ‘, u w,hl *) 
has a positive lower bound uniformly. Combing (14.361) . (14.371) . (14.391) . (14.401) 
and (13351) yields (TP2D . □ 

Remark 4.1. fli in (13.71) and (13.91) can be different from that in (13.61) and 
(13.81) . which also ensure that the corresponding estimates in Theorem 4.2- The¬ 
orem 4.4 still hold. 


Remark 4.2. By dropping the steps computing u 


Hi ' 


, e h *,u w ’ hi *,\ w ’ki¬ 


rn Scheme 3.2, and by replacing X w = 

\w _ a(u w ,u w ) i \w,hi _ a(u w ’ hi ,u w h "' 
b(u w ,u w ) a(u w ’ h i ,u u 




l and \ w ’ hi * = ^ 


*) 


with 


77 -y, respectively, we are able to establish 
multilevel discretizations based on local defect-correction for symmetric Eigen¬ 
value Problems. Hence the corresponding estimates in Theorem 4.1-Theorem 
4.4 still hold. 


Remark 4.3. By referring to [ 8 ], we can establish the parallel version 
of Scheme 3.1 and Scheme 3.2 and have the corresponding error estimates in 
Theorem 4.3-Theorem 4.4 apparently. 


5 Numerical experiments 

5.1 Computational method for (A * H ,u* H ) (see | I24L l25 j) 

Assume that (Ah , uh) is obtained from Scheme 3.1 or Step 1 of Scheme 3.2, then 
Ah = Ah, and from |24 [ 125 ) we can obtain u* H by using the following approach 
such that \b(un,u* H )\ has a positive lower bound uniformly with respect to H. 
Let mo be the algebraic multiplicity of Ah and l be ascent of Ah- 
Let ujf be the orthogonal projection of uh to N((jk -T h) 1 ), and u* H = 

u h/\\ u h Ho- When uh € N((-p - Th) 1 ) it is clear that u* H = uh- When 

uh A5((—h— T^Y), to find u * H , first we seek a basis {4>i}T° °f —2 
and solve the following equations 


m 0 


^ ^ {$i ? 4*j ) — ( U H •> 4*j ) 5 3 — 15 2 , • ’ * , TTIq , 

(5.1) 

2=1 


mo 


U H = ^ ^ 0^2 02 ; 

(5.2) 

2=1 


U H = u h/\\ u h\W 

(5.3) 


Obviously, u* H satisfies 

(ti Tff) l u* H = 0, |Kllo = l. 

a h 


16 








































Thus, to find u* H which satisfying (15.21) and (15.31) in V£(fi), the key is to 
seek a basis of N((j7 — ThY)- 

When Z = 1, it is actually to solve the following equations to obtain a basis 
in the solution space. 

u [ h ] g vg(n), 

a(v, u(jp) - X H b(v, u$) = 0, Vv <E V$(f2). (5.4) 

(When Xh is a simple eigenvalue, l = 1 and —Tj}) z ) is a one-dimensional 
space spanned by the eigenfunction u^.) 

When l > 1, how to seek a basis {</>}of N((-^ - T^) 1 ) efficiently is an 

important issue of linear algebra. 

5.2 Numerical Examples 

Consider the convection-diffusion equation 

— Au + b ■ Vu = Xu, in 12, u = 0, on <9f2, (5-5) 

where Q = (—1, l ) 2 \ {[0,1] x [—1,0]} or 12 = (—1,1) 2 \{{0} x [—1,0]}. The 
first eigenfunctions of both problems have the singularities at the origin. The 
exact eigenvalues, which are unknown, are thereby replaced by approximate 
eigenvalues with high accuracy. For the problem with b = (1,1) T , b = (0,3) T 
and b = ( 0 , 10 ) T on 12 = (— 1 , 1) 2 \ {[ 0 , 1 ] x [— 1 , 0 ]}, we take the approxi¬ 
mate first eigenvalue as Ai ss 11.8897, Ai « 10.1397 and Ai « 34.6397, re¬ 
spectively. For the problem with b = (1,1) T , b = (0, 3) T and b = (0,10) T 
on 12 = (— 1 , 1 ) 2 \{[ 0 , 1 ] x { 0 }}, we take the approximate first eigenvalue as 
Ai ~ 10.621,Ai ~ 8.871 and Ai « 33.371, respectively. We will report some nu¬ 
merical experiments by using linear finite elements on uniform triangle meshes. 
In our numerical experiments, we use Scheme 3.2 to solve the problem such that 
= ({jr, Tp -) 2 \ {[0, i-] x [=r, 0]} for L-shaped domain, 12; = (=£ , ) 2 \{{0} x 

[—^r, 0 ]} for slit domain, i = 1 , 2 , ••• , 6 , and locally fine grids have the same 
degree of freedom as that of globally mesoscopic grid (see Tables 1-6). 

In our experiments, according to the assumptions of Theorem 4.4, we ap¬ 
proximately take 71 = 72 = 1/2, 2/3 and s = S 2 = 1/2, 2/3 so that (14.351) holds 
for slit domain and L-shaped domain, respectively. We use MATLAB 2011b 
under the package of Chen (see 1 ) to solve the problems, and the numeri¬ 
cal results are shown in Tables 1-6. From Tables 1-4 we can see that without 
increasing degree of freedom on locally fine grids, the first local defect correc¬ 
tion can largely improve the accuracy of the eigenvalues, and the local defect 
corrections that follows can gradually improve the accuracy of the eigenvalues 
by overcoming the singularity at the origin. But Tables 5-6 also indicate that 
Scheme 3.2 is not valid for the problems with b = (0,10) T . Concerning this 
point, the figures of eigenfunction and its adjoint pair (see Fig. 5.1-5.2) shows 
their functional-value abrupt changes mainly center on boundary layer, which 
may lead to the invalidity of Scheme 3.2. for the case b = (0,10) T , we adopt 
the parallel version of Scheme 3.2 to make local defect-corrections on boundary 
layers with functional-value abrupt changes (see also Fig. 5.1-5.2). 

Specifically speaking, for the L-shaped domain, we find that it’s better 
to make local defect-corrections near the origin on slightly small area fl} = 
(^irr, x [— 2 TTT 7 0 ]} f° r both eigenfunction and its adjoint eigenfunc¬ 

tion; as for the other local defect-correction areas, we set as fl 2 = (^v/t; yr/r) x 
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(1 - 1) for the eigenfunction, fif = (-1 - \ + jr) x (—1, —1 + d.) for 

the adjoint eigenfunction, respectively; the related numerical results is given in 
Table 7. Here we set 

DOF w « | x DOF n i « | x DOF Q 2 «4x (i = 1,2,---)- 

For the slit domain, we set as the local defect-correction area 12} = ^r) 2 \{{0} x 

[— 0 ]} for both eigenfunction and its adjoint eigenfunction, f 2 ? = ^r) x 

(1 — i-, 1 ) for the eigenfunction, 12 ? = (| — dj, A + dj.) x (— 1 , —1 + i.) and 
12 } = (—5 — — \ + 5 J-) x (— 1 , —1 + d,) for the adjoint eigenfunction, respec¬ 

tively; the related numerical results is given in Table 8 . Here we set 

DOF w = DOF n i = DOF n 2 « 2 x DOF n 3 = 2 x DOF n * (* = 1,2,---). 

Table 7.12 and Table 7.16 in [5] show that, using the adaptive homotopy 
method to solve the L-shaped domain problem with b = (10,0) T , the ap¬ 
proximate eigenvalue can have 4-5 significant digits with DOF = 154994 and 
124469, thus the adaptive homotopy method is efficient. However, by using 
our algorithm, the approximate eigenvalue can have 6 significant digits with 
DOFh = 12033 (see Table 7), which also indicates our algorithm is efficient. 


Table 1: 

12 = (-l,l) 2 \{[0,l]x 

[-1,0]}, b 

= (0,3) T . 

DOF h 

DOF w 

A H 

X w 


705 

2945 

11.94916 

11.91247 

11.89949 

2945 

12033 

11.91250 

11.89859 

11.89343 

12033 

195585 

11.89859 

11.89109 

11.89028 

^ W,h2 


yw,h4 

X w,h$ 


11.89466 

11.89275 

- 

- 

- 

11.89146 

11.89068 

11.89037 

~ 

- 

11.88996 

11.88983 

11.88978 

- 

- 

Table 2: 

12 = (-l,l) 2 \{[0,l]x 

[-1,0]}, b 

= (Mf. 

DOF h 

DOF w 

A H 

x w 

yw,h\ 

705 

2945 

10.21836 

10.16730 

10.15332 

2945 

12033 

10.16730 

10.14979 

10.14438 

12033 

195585 

10.14979 

10.14117 

10.14034 

\ w ,h2 

\w,h 3 

yw,h4 

\w,K 

\w,h 6 

10.14836 

10.14643 

- 

- 

- 

10.14238 

10.14160 

10.14129 

- 

- 

10.14002 

10.13989 

10.13984 

- 

- 
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Table 3: 

fi = (-l,l) 2 \{{0}x [-1,0]}, b 

= (0,3) T . 

DOF h 

DOF w 

A H A™ 



945 

3937 

10.79397 10.70640 

10.66393 

3937 

16065 

10.70630 10.66356 

10.64247 

16065 

64897 

10.66353 10.64237 

10.63186 

yw,h 2 

\w,h 3 

^w,h 4 yw,h*, 


yv,h e 

10.64333 

10.63315 

- 


- 

10.63208 

10.62691 

10.62434 - 


- 

10.62664 

10.62404 

10.62274 10.62209 

~ 

Table 4: 

0 = (-1,1) 2 \{{0}x [-1,0]}, b - 

= (M) T . 

DOF h 

DOF w 

A H A™ 

X w 

,hi 

945 

3937 

9.06244 8.96099 

8.91741 

3937 

16065 

8.96090 8.91470 

8.89333 

16065 

64897 

8.91468 8.89266 

8.88207 

^w,h 2 

\w,h 3 

^w,h 4 ^W,/l5 

X w 

,h*Q 

8.89663 

8.88641 

- 

- 


8.88289 

8.87772 

8.87514 - 

- 


8.87684 

8.87424 

8.87294 8.87229 

~ 


Table 5: 

^ = (-1,1) 2 \{[0,1] x [1,0]}, 6 

'=(0,10) T . 

DOF h 

DOF w 

A H X w 


^w,h i 

705 

2945 

34.58756 34.63484 

34.61999 

2945 

12033 

34.63473 34.64168 

34.63605 

12033 

195585 

34.64167 34.64066 

34.63982 

^w,h 2 

\w,h 3 

yw,hi \w,h 5 


\w,h 6 

34.61632 

34.61460 

- 


- 

34.63437 

34.63363 

34.63333 - 


- 

34.63952 

34.63939 

34.63934 - 


- 


Table 6: 

fi = (-l,l) 2 \{{0}x 

[-1,0]}, b 

= (0,10) T . 

DOFh 

DOF w 

A H 

X w 

^w,hi 

945 

3937 

33.43287 

33.42950 

33.38763 

3937 

16065 

33.42885 

33.40686 

33.38587 

16065 

64897 

33.40671 

33.39070 

33.38021 

)W,h 2 

\w,h 3 

^w,h 4 

\w,h 5 

\w,h 6 

33.36894 

33.35923 

- 

- 

- 

33.37595 

33.37090 

33.36836 

- 

~ 

33.37510 

33.37253 

33.37124 

33.3706 

~ 
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Table 7: 

fi = (-l,l) 2 \{[0,l]x 

[-1,0]}, b 

= (0,10) T . 

DOF h 

DOF w 

A H 

x w 

^w,hi 

705 

2945 

34.58756 

34.63484 

34.64245 

2945 

12033 

34.63474 

34.64169 

34.64162 

12033 

195585 

34.64166 

34.64067 

34.64017 

^W,h,2 



\^w,h 5 

}w,h 6 

34.64102 

34.63949 

- 

- 

- 

34.64049 

34.63980 

34.63951 

~ 

~ 

34.63990 

34.63978 

34.63973 

- 

- 

Table 8: 

fi = (-l,l) 2 \{{0}x[ 

-1,0]}, b = 

= (0,10) T . 

DOF h 

DOF w 

A H 

A w 

^w,h 1 

945 

3937 

33.43287 

33.42950 

33.40090 

3937 

16065 

33.42885 

33.40686 

33.38916 

16065 

64897 

33.40671 

33.39070 

33.38103 

}w,h 2 

yw,h 3 

^w,h 4 

\w,h 5 

yw,h 6 

33.38429 

33.37474 

- 

- 

- 

33.37976 

33.37475 

33.37221 

- 

- 

33.37605 

33.37349 

33.37220 

33.37155 

- 



Fie. 5.1. Eigenfuntion and its adjoint pair with b = (0, 10) T on = 
( I 51 ) 2 \{[ 0 , 1] x [—1.0]} 



Fig. 5.2. Eigenfuntion and its adjoint pair with b = (0, 10) T on 12 = 
(-1,1) 2 \{{0} x [-1,0]} 
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